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SDSS J1130+0058 an X-shaped Radio Source With 
Double-Peaked Low-Ionization Emission Lines: A binary 
Black Hole System? 
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ABSTRACT 

In this paper we study the object SDSS J1130+0058 which is the only AGN known 
to have both double-peaked low-ionization broad emission lines, and also X-shaped 
radio structures. Emission from an accretion disk can reproduce the double-peaked 
line profile of broad Ha, but not the radio structure. Under the accretion disk model, 
the period of the inner emission line region is about 230 years. Using a new method to 
subtract the stellar component from the data of the SDSS DR4, we obtain an internal 
reddening factor which is less than previously found. The implied smaller amount of 
dust disfavors the backflow model for the X-shaped radio structure. The presence of a 
Binary Black Hole (BBH) system is the most natural way to explain both the optical 
and radio properties of this AGN. Under the assumption of the BBH model, we can 
estimate the BBH system has a separation of less than 0.04 pc with a period less than 
59 years, this may pose some problem to the BLRs sizes, still we conclude that the 
BBH model is favored on the basis of the present limited information. 

Key words: Galaxies: Active - galaxies muclei - galaxies:jets - accretion, accretion 
disks - black hole physics 



1 INTRODUCTION 

The detailed structure of the inner parts of Active Galactic 
Nuclei (AGN): the black hole, accretion disk and broad emis- 
sion line regions (BLRs) , cannot be resolved by direct obser- 
vations, with the possible exception of indirect evidence from 
a dusty torus in the nucleus of NGC4261. The use of rever- 
beration mapping (Peterson 1993; Netzer & Peterson 1997; 
Blandford & Mckee 1982) requires the assumption of a rough 
geometric structure for the BLRs in order to explain the cor- 
relation between the variations of broad emission lines and 
those of continuum. Moreover, the problem of spatially re- 
solving these regions, will not be solved in the foreseeable 
future. Up until now, the best approach is the detailed study 
of broad hue profile shapes and asymmetries from high S /N 
spectra (Sulentic et al. 2000; Marziani, Dultzin-Hacyan & 
Sulentic 2006). 

There are at least three special types of AGN where the 
observed phenomenology has suggested the presence of a bi- 
nary black hole system (BBH) in the center: (1) blazars with 
periodical brightness variations, (2) AGN with X-shaped 
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radio structures (X-shaped AGN), and (3) double-peaked 
broad low-ionization line (dbp) emitters. 

In the prototypical case of the periodically varying 
blazar OJ 287 (Sillanpaa et al. 1996a, 1996b; Valtonen et 
al. 2006a) the bending of the VLBI (Very Long Baseline In- 
terferometry) jet was first reported by Vicente et al. (1996). 
A very small change in the orientation of the jet is needed to 
change the Doppler boosting dramatically, thereby produc- 
ing long-term periodic brightness modulations. Sillanpaa et 
al. (1988) and Lehto & Valtonen (1996) modeled the peri- 
odic outbursts by relating them to tidally induced mass flows 
from the accretion disks in a BBH system with an orbital 
period of 8.9 years (in the rest frame of OJ 287). The predic- 
tions of the model have been confirmed quite spectacularly 
(Valtonen et al. 2006a; Valtonen et al. 2006b). 

The first and most famous dbp emitters are NGC 
1097 (Storchi-Bergmann, Nemmen, et al. 2003; Storchi- 
Bergmann, Eracleous et al. 1997; Storchi-Bergmann, Era- 
cleous & Halpern 1995; Storchi-Bergmann, Baldwin et al. 
1993), Arpl02B (Chen et al. 1989, 1997; Chen & Halpern 
1989; Halpern et al. 1996; Antonucci et al. 1996; Sulentic et 
al. 1990) and 3C390.3 (Shapovalova et al. 2001; Gilbert et 
al. 1999). The first natural explanation for the appearance of 
double peaked emission lines was also a BBH model (Begel- 
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man et al. 1980; Gaskell 1983). Another model of bipolar 
outflow was proposed by Zheng et al. (1990). The most suc- 
cessful model involved lines originating from an accretion 
disk to explain not only the appearance, but also the ob- 
served variations of doublc-pcakcd emission lines (Chen et 
al. 1989; Chen & Halpern 1989). The simple accretion disk 
model has been modified to better fit the observations to 
an elliptical disk by Eracleous et al. (1995), to a warped 
accretion disk by Bachev (1999) and Hartnoll & Blackman 
(2000), and to a circular accretion disk plus spiral arms by 
Chakrabarti & Wiita (1994), Hartnoll & Blackman (2002) 
and Karas, Martocchia & Subr (2001). An important reason 
to give up the initial BBH model was the unreasonable BH 
masses derived from the variations of the broad emission 
lines for three dbp emitters (Eracleous et al. 1997). In the 
case of the bipolar outflow model, the observed variations 
could not be explained. 

Another special kind of AGN, X-shapod AGN, are clas- 
sified according to their extended radio morphology. One of 
the first explanations for the distortion of radio structures, 
particularly for X-shaped and Z-shaped AGN, was backflow 
material dominating the bridge regions (Leahy & Williams 
1984; Worrall et al. 1995). More recently, other models have 
also been developed to explain the distortions, such as the 
precession of jot axis in a BBIf system. It should be no- 
ticed that in the case of OJ287 (the first confirmed peri- 
odically varying blazar), the bending of the VLBI jet was 
first reported by Vicente et al. (1996). Objects with bend- 
ing, misalignment, and wiggling of extragalactic radio jets 
(often associated with knots superluminally moving along 
different-scale curved trajectory), have been interpreted in 
terms of helical structures of the jets. This structure is likely 
caused by the precession of the jet in a binary black hole 
system (both kinds of models are discussed in e.g.: Dennett- 
Thorpe et al. 2002; Merritt & Ekers 2002; Liu 2004; Caproni 
& Abraham 2004a, 2004b; Ostorcro, Villata & Raitcri 2004). 
Up to the present, there arc more than ten X-sliapcd AGN 
reported in the literature. 

It seems only natural to invoke a BBH system to explain 
the presence of both dbp lines and X-shaped jets in an AGN. 
In this paper, we report the only such object so far known: 
SDSS J113021.4H-005823 (hereafter, SDSS J1130+0058, 
found in SDSS DR4 (York et al. 2000; Strauss et al. 2002; 
Abazajian et al. 2004; Berk et al. 2001; Scott et al. 2004)). 
The next section presents the observational data. In Sec- 
tion 3 we present the discussion and conclusions. In this pa- 
per, the cosmological parameters Ho = 70km • s~^Mpc~'^, 
Qa = 0.7 and fim = 0.3 have been adopted. 



2 OBSERVED RESULTS 
2.1 Results at Optical Band 

SDSS Jl 130+0058 has been studied by Wang et al. (2003) 
with another target name (4C +01.30), as an obscured 
quasar. Thus, in this section, we briefiy discuss the observed 
results at optical band. The apparent optical Petroson mag- 
nitudes of SDSS J1130+0058 {z = 0.13) are 18.43, 17.13, 
16.19, 15.72, 15.49 at u, g, r, i and z bands, respectively. The 
galactic reddening corrected (with E(B-V) = 0.022) spec- 
trum of SDSS Jl 130+0058 with total exposure time of 3600 
seconds is shown in Figure 1. 



In order to measure the line parameters more accu- 
rately, the component due to the host galaxy must be re- 
moved from the observed spectrum. We believe that the 
most convenient way to subtract this component is the PCA 
(Principle Component analysis) method described by Li et 
al. (2005) and Hao et al. (2005), using the eigenspectra from 
pure absorption galaxies from SDSS or the eigenspectra from 
stars in STELIB (Le Borgne et al. 2003). Here, we used the 
method from Hao et al. (2005). The eigenspectra are calcu- 
lated by KL (Karlmnen-Locvc) transform for about 926 pure 
absorption galEixies selected from SDSS DR2. Then, the first 
eight eigenspectra and the spectrum of an A star are used 
to fit the absorption properties of the observed spectrum. A 
power law is used to fit the contirmum from the rmclei. The 
line spectrum of SDSS J 1130+0058, after the subtraction of 
the stellar component and power law continuum, are show in 
Figure 2. The stellar components and the power law contin- 
uum are also shown in the figure. The best fitted results for 
absorption lines, such as the CaIIA3934, an6i3974A doublet, 
indicate that the subtraction of stellar components is reli- 
able. As discussed in Wang et al. (2003), there is apparent 
broad Ha emission line and weak broad H/3 emission line in 
line spectrum. 

After the subtraction of the stellar light and contin- 
uum, the emission line parameters are measured. The line 
parameters were listed in Wang et al. (2003), however, 
new eigenspectra and new observed spectrum from SDSS 
DR4 are used here, thus, we re-measured the emission 
lines and list our new line parameters in Table 1. We fit 
each narrow emission line by a narrow gaussian function. 
For [OIII]A4959,5007A doublet, there are other two gaus- 
sian functions for each line, because of the asymmetry of 
[OIII]A4959, 50071. Due to the weakness of broad H/3, we 
use the broad component of Ha to fit broad H/J. First, we 
obtain the broad component of Ha after the subtraction of 
narrow Ha, [NII]A6548, 65831 and [SII]A6716, 673ll within 
the wavelength range from 6400l to 68OOA. In order to get 
more accurate narrow components near Ha, there are two 
broad gaussian functions to fit the broad Ha. Then the nar- 
row components near Ha according to the best fitted results 
arc subtracted from the observed line spectrum. We have as- 
sumed that broad H/3 and Ha have the same line profiles, 
thus we can estimate the fiux ratio of broad Ha to H/3. The 
best fitted results for Ha are shown in Figure 3. Figure 4 
shows the best fitted results for H/3 by scaling of the ob- 
served broad Ha. 

The object SDSS J1130+0058, has been considered as 
a dbp emitter in the sample of dbp emitters by Strateva et 
al. (2003). In our new sample of dbp emitters (Zhang et al. 
in preparation) from SDSS DR4, this object also is a dbp 
emitter. For preliminary definition of criteria see also Zhang, 
Dultzin-Hacyan & Wang (2007). 

2.2 Results at Radio Band 

The radio image at 20cm for SDSS J1130+0058 can be 
extracted from Faint Images of the Radio Sky at Twenty 
centimeter Survey (FIRST, Becker, White & Helfand 1995 
), which is shown in Figure 5. The apparent X-shaped ra- 
dio morphology can be seen. There are four regions which 
strongly emit radio power as listed in Table 2. As discussed 
in Wang et al. (2003), the smallest offset of the photometric 
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center from the four radio regions is about 4.5 arcseconds 
which is larger than the position uncertainty of SDSS and 
FIRST, which indicates that the central component may be 
a component from the jet rather than the radio core. The 
total integral radio power at 20cm for SDSS J1130+0058 is 
about 1025'*^W-Hz-\ 

Furthermore, one fundamental parameter of a radio jet, 
the age of the jet, is important to understand the physical 
processes in the forming the jet. There are many papers that 
discuss the age of a jet, including the dynamical and spectral 
age. Alexander & Leahy (1987) firstly proposed the method 
to measure the spectral age of a radio jet (the same methods 
can be found in Leahy ct al. 1989; Liu ct al. 1992; Params et 
al. 1999). Spectral age depends on the break frequency Hbr, 
which can be calculated by the spectral index at different 
frequencies in various regions. The more accurate dynamical 
age is more difficult to calculate and depends on the model 
of dynamical expansion (Kaiser & Alexander 1997; Kaiser 
2000; Kaiser ct al. 1997; Parma et al. 1999; Machalski et al. 
2006; Tudose et al, 2006; Kino & Kawakatu 2005). The hmit 
on the radio data is not enough to calculate the accurate age. 
However, we can determine which one is older. Machalski et 
al. (2006) have found there is a strong linear correlation 
between spectral age and dynamical age (dynamical age is 
about two times the spectral age). The spectral age ts can 
simply depend on the distance between the lobe edge and 
the core x: ts oc x, considering the relation between x and 
Hijr shown in Parma ct al. (1999). Tims as a simple empirical 
result, we can say that the double-jet in the N-S direction is 
older than the one in E-W direction, because of the longer 
projected distance in N-S direction. 



3 DISCUSSIONS AND CONCLUSIONS 
3.1 An obscured AGN? 

The internal reddening factor can be determined by the 
value of the Balmer decrement. The flux ratios of the Balmer 
emission lines are about 3.78 and 6.52 for the narrow and 
broad components respectively. It should be pointed out, 
however, that the intrinsic Balmer decrement is not a defini- 
tive value for AGN, because case B recombination is not 
necessarily valid (particularly for the BLR). The value of 3.1 
can be expected by Case B recombination. From a sample of 
quasars with smaller contribution from the host galaxy, the 
mean value is about 3.5 (Greene & Ho, 2005b). If we assume 
the intrinsic Balmer decrement as 3.1, the internal redden- 
ing factors are about E(B-V)=0.64 for the broad Balmer 
emission lines, and E(B-V) = 0.17 for the narrow Balmer 
emission lines. The Balmer decrement is not as large as the 
one determined by Wang et al. (2003). 

In order to verify if the internal reddening factor leads 
to luminosities which agree with well established scaling 
relations, we have used the tight correlation between con- 
tinuum luminosity at 5100A and the luminosity of Ha 
found by Greene & Ho (2005b). The observed luminos- 
ity of Ha including the narrow component of Ha is about 
1.16 X 10''^ erg • s"^ ( 9.08 x 10*^erg-s"^ for broad compo- 
nent and 2.47 x lO^^erg • s^^ for narrow component). The 
continuum luminosity at 5100A from the nucleus can be 
measured from the power law continuum, L = = 1.09 x 



10*'^ erg ■ s^^. The total continuum luminosity at 5100A from 
the observed spectrum is approximately 4.79 x lO'^'^erg • s~^. 
After the internal reddening correction with E{B — V) = 
0.64, the intrinsic continuum luminosity at 5100A from the 
nucleus is about 8.03 x lO^^crg ■ s^^, and the intrinsic lumi- 
nosity Lhc, is about 4.01 x 10''^ erg • s~^ (the luminosity of 
narrow Ha has been corrected by the factor E(B-V)=0.17). 
The relation between the corrected values are consistent 
with the expectation from the correlation found by Greene 
& Ho (2005b): 

= 5.25 X 10«( JM^)"-erg • s^ (1) 

This consistency confirms that the method of subtract- 
ing stellar components is valid and the calculated internal 
reddening factor is a reasonable estimate. SDSS J1130-I-0058 
is indeed an obscured AGN with intermediate internal red- 
dening. The smaller value for internal reddening found in 
this paper (as compared to Wang et al. 2003) implies a 
smaller amount of internal dust, and thus a lower density 
for the internal regions. However, the backflow model re- 
quires a higher gas pressure in the direction of the galaxy, 
as shown in other X-shaped radio sources there are nearly 
no detectable broad emission lines because of the higher in- 
ternal reddening. This does not favor the backflow model for 
the X-shaped radio structure in SDSS J1130-I-0058. 



3.2 BH Masses and Accretion Rate 

There are mainly two methods to estimated the central BH 
masses: One is from the correlation between the masses of 
the bulge and the central Black Hole following the strong 
correlation between BH masses and stellar velocity disper- 
sion (Tremaine et al. 2002; Ferrarese & Merritt 2001; Geb- 
hardt. Bender ct al. 2000) (or line width of narrow emission 
lines, Greene & Ho 2005a; Nelson & Whittle 1995, 1996; 
Nelson 2000), and the correlation between absolute magni- 
tude of the bulge and BH masses (Kormendy & Gebhardt 
2001) etc. The other one is from the assumption of virializa- 
tion (Kaspi ct al. 2000; Peterson & Wandol 1999; Peterson 
et al. 2004; Onken et al. 2004; Sulentic et al. 2006), based 
on the line width of broad emission lines and the size of 
broad emission line regions estimated from continuum lumi- 
nosity (Kaspi et al. 2000; 2005). As shown in Kazantzidis 
ct al. (2005), the central masses of merging galaxies should 
also obey the strong correlation between BH mgisses and 
stellar velocity dispersions, if they undergo a dissipational 
process with star formation mergers. Moreover, according to 
the hierarchical merging model, wc think the BBH system is 
common in the center of AGN. Even in the sample used to 
obtain the strong relation Mbh — o (Tremaine et al. 2002; 
Ferrarese & Merritt 2001; Gebhardt, Bender et al. 2000), we 
cannot discard the possibility that the centers of the objects 
have BBH system. Thus, whether there is a BBH system in 
the center of SDSS J1130+0058 or not, we estimate the BH 
masses from stellar velocity dispersions. 

How to measure accurately stellar velocity dispersions 
is an open question, because of the known problem about 
the template mismatch. A commonly used method is to 
select spectra of several kinds of stars (commonly, G and 
K) as templates, and then broaden the templates by the 
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same velocity to fit stellar features, leaving the contribu- 
tions from different kinds of stars as free parameters (Rix 
& Whit 1992). However, more information about stars in- 
cluded by the templates should lead to more accurate mea- 
surement of stellar features. According to the above men- 
tioned fitting method, we selected a new template rather 
than several spectra of G or K stars as templates. The 
method of Principle Component Analysis (PCA) provides 
a better way to constrict more favorable information from 
a scries of spectra of stars into several eigcnspcctra. Thus, 
we apply PCA method for 255 spectra of different kinds 
of stars in STELIB. Selecting the first several eigenspectra 
and a three-order polynomial function for the background 
as templates, the value of the stellar velocity dispersion can 
be measured by means of the method of minimum fit. 
Figure 6 shows the fitted results for the stellar features near 
MgIbA5175A by means of the eigenspectra with broadening 
velocity a = 142.6 ± 22.6km • s"^ The error in the stellar 
velocity dispersion is determined from the different numbers 
of required eigenspectra. Here, 4 to 8 eigenspectra are used 
to fit the stellax features. The central BH masses of SDSS 
J1130-I-0058 can be estimated as Mbh = S.SStt'tj x lO^'M® 
by (Tremaine et al. 2002; Ferrarese & Merritt 2001; Geb- 
hardt. Bender et al. 2000): 



Me 



-|^q8.13±0.06^ 



200km • s 



n4.02±0.32, , 



(2) 



Moreover, we should select another parameter to estimate 
the central BH masses. The convenient way is to select the 
line width of narrow emission lines to trace the stellar veloc- 
ity dispersion in the bulge. From the line widths of narrow 
emission lines listed in Table 1, the minimum value of line 
width {aiine) is 103.3km • s~^ for H(5 and the maximum value 
is 175.6km • s^^ for [OIIJA3727A. Thus, we can estimated the 
BH masses as 9.47 x 10® - 8.01 x lO'^Mg by the relation of 
Mbh — cr, which is consistent with the value from stellar 
velocity dispersion and confirms the strong correlation be- 
tween stellar velocity dispersion and line width of narrow 
emission lines. 

The virial BH masses are based on the empirical re- 
lation between the size of BLRs and continuum luminos- 
ity. However, the broad emission lines of SDSS J1130+0058 
have special double-peaked profiles. Thus, either under the 
accretion disk model or under the BBH model, the size of 
BLRs of the object should not obey the empirical relation 

For accretion rate, we can determine the dimensionless 
accretion rate based on the bolometric luminosity and Ed- 
dington luminosity to trace the accretion rate. The bolo- 
metric luminosity for AGN with normal spectral energy 
distribution, such as normal big blue bump, can be de- 
termined by the continuum luminosity Lhoi ~ 9 x L^^^^^. 
The Eddington luminosity can be determined by LEdd 



1. 38 X 10^* Mbh /Mc3' 



Thus, the dimensionless accre- 



tion rate can be determined as rn ~ Lhoi / Lsdd ^ 0.06 — 0.36, 
which is larger than the critical accretion rate for ADAF 
(Mahadevan & Quataert 1997; Narayan et al. 1995, 1996; 
Mahadevan 1997). 



3.3 Accretion Disk Model or BBH model for 
double-peaked Ha? 

First, we check whether the BBH model can be valid for 
SDSS J1130-I-0058. From the best fitted results shown in 
Figure 3, double-gaussian functions can completely repro- 
duce the observed double-peaked broad Ha. If the BBH 
model is valid, the observed radial velocity of each peak 
of each broad component of Ha represents the rotating ve- 
locity in the line of sight, ui ~ 1065km ■ s~^ and V2 ~ 
— 2354km • s"'^. Here, the minus means the component is 
moving toward us. Then the mass ratio of the two central 
BH masses can be roughly determined by Mbhi/Mbh2 ~ 
vi/v\ ~ 2.2, if the separation between the central two black 
holes is much larger than the distance between BLRs and 
each corresponding black hole. Here, Mbhi represents the 
black hole with larger BH mass. However, whether the re- 
lation Mbh — c7 is valid for merging galaxies depends on 
whether the merger is dissipative or colisionless (Kazantzidis 
et al. 2005). Thus, we determine the BH masses using an- 
other method, before considering the total BH masses from 
stellar velocity dispersion. We can estimate the BH masses of 
central two black holes under the correlation Mbh — L 

5100A 

found by Peterson et al. (2004), Mbhi ~ 6.5 ± 1.1 x lO^M© 
and Mbh2 ~ 3.1 ± 0.9 x IO^Mq, which is consistent with 
the mass ratio from the ratio of rotating velocity. The con- 
tinuum luminosity of each black hole -|- BLRs system can be 
estimated from the luminosity of each broad component of 
Ha by equation 1, because of the small contributions from 
narrow components of Ha. Under the binary black hole sys- 
tem model, there is another parameter of projected angle of 
the orbital velocity to determine the structure of each BLR 
except the inclination angle of the orbit, thus we select the 
correlation between Mbh — L - rather than eq 2. The 
total BH masses, 9.6 ± 2.0 x IQ'M.q, determined from each 
broad Ha is similar to BH masses 3.83^2 37 x 

es- 
timated from stellar velocity dispersion, which is consistent 
with the predicted result for merging galaxies by Kazantzidis 
et al. (2005). 

Under the assumption of the BBH model, we accept the 

BH masses estimated from each broad Ha for each black 
hole. We can roughly estimate the separation of the central 
two black holes by: 



,2 _ G X {Mbhi + Mbh2) 



wi = O X 



Ainu. 



(3a) 



Mbhi + Mbh2 



X r X sin(i) sin(0) (3b) 



where r is the separation between the central two black 
holes, i is the inclination angle of the rotating plane, 4> is 
the projected angle of the orbital velocity, Q is the angular 
speed, Vi is the observed rotating velocity in the line of sight 
of the clouds in BLRs including Mbhi- From the equations 
above, we can determine the upper limit value of the sep- 
aration r: r < 40.8 light — days, and the upper limit value 
of orbital period P: P < 59 years. However, the size of each 
BLR can be determined by (Wang & Zhang 2003): 



logRBLRs =(0.51 ± 0.02) log( 



lO-i^erg • s-i ' 
-I- (2.16 ± 0.02)light - days 



(4) 



The sizes of BLRs are about 45.1 light-days and 79.4 
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light-days for each BLR, which turns out larger than the 
maximum value of separation. In order to keep the com- 
pleteness of each BLR, the sizes of the BLRs should not be 
larger than the separation of the two black holes. This poses 
some problem to the BBH model for SDSS J1130+0058. It 
has to be stressed, however, that the standard way to deter- 
mine BLR sizes, may not be applicable for a BBH system. 

Another model, the accretion disk model has been ap- 
plied to explain the double-peaked Ha. Here, we select ellip- 
tical accretion disk model (Eraclcous ct al. 1995), because 
of the number of free parameters in the model which is less 
than those required by the circular accretion disk plus spi- 
ral arms model. The best fitted results with ~ 0.91 are 
shown in Figure 7. The inner radius is about 713 ± 107Rg, 
outer radius is about 6327 ± 939Rg, eccentricity is about 
0.22±0.01, the inclination angle of the disk is about 39±3°. 
The slope of line emissivity is about e oc The lo- 

cal broadening velocity is about 539 ± 56km ■ s~^. According 
to the eccentricity, the period of the inner emission regions 
can be estimated as Ppre ~ 10.4 ^-^^^^^1.5 ^^61^3^ ^ 230years, 
where Ms is the central BH masses in unit of lO^M©, rs is 
the radius in unit of lO'^Rc;. 

More information, particularly about the variations 
of double-peaked broad Ha to determine the time scale 
and other characteristics of the variations of the lines is 
needed to verify whether the BBH model is valid for SDSS 
J1130+0058. From the theoretical side, a detailed model of 
the combined accretion disk emission from two accretion 
disks is yet to be developed. 



3.4 Final Considerations 

The only X-shaped radio source with apparent double- 
peaked broad emission lines, SDSS J1130-f0058, is indeed 
an internally obscured AGN with standard contirmum lumi- 
nosity of a Seyfert 1 galaxy. The presence of internal dust 
may provide a convenient background for the backflow of 
material to produce X-shaped radio structures. But in this 
paper we found a smaller value of the reddening (and thus 
of the dust contents) than previously reported. We are not 
in the position to quantitatively support nor reject the back- 
flow model, however a natural alternative is the BBH model. 
Objects with X-shaped extragalactic radio jets may be the 
signature of these binary systems. 

A BBH system can also account for the observed dbp 
Balmer emission lines. Except for the problem with the sizes 
of the BLRs, the separation between the two black holes un- 
der the assumption of BBH model, is consistent with the the- 
oretical expectation by Merritt & Ekers (2002), who found 
that the separation of the two black holes should be about 
0.01 pc to 1 pc in a time shorter than 1 Gyr, if the co- 
alescence rate of binary black holes is comparable to the 
galaxy merger rate. Under the BBH model, the separation 
of central black holes for SDSS ,11130-1-0058 is about less 
than 0.04 pc. From the results of a single observation, we 
cannot rule out the BBH model. Moreover, the periods of 
BBH model and accretion disk model are very different, one 
is about 230 years and the other less than 59 years, so that 
long term observations are necessary to understand SDSS 
Jl 130-1-0058 better. Furthermore, the BBH model and ac- 
cretion disk model will predict different variation patterns 



of broad Ha. The BBH model leads to radial velocity of 
each peak varying in period. On the other hand, the accre- 
tion disk model cannot lead to the same result. Long term 
monitoring should provide more information to discriminate 
which model for SDSS J1130+0058 is favored. 

There is another possibility within the BBH system sce- 
nario to explain the double-peaked broad emission lines and 
the X-shaped radio structures. The double-peaked broad 
emission lines are coming from the accretion disk of only one 
of the two central black holes. There is NO BLR around the 
other black hole. Otherwise, there would be three or more 
components of broad Ha. The X-shaped radio-structures are 
produced by the precession of the jet. The main effect of the 
existence of the second black hole is to produce the preces- 
sion of the radio jet with respect to the accretion disk. In this 
case, the predicted behavior is that the center wavelength 
of broad Ha should vary in period, besides the variation of 
line profiles. Long term observations also can provide infor- 
mation to distinguish this model from the other two models. 
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Table 1. Line Parameters of SDSS J1130+0058 



Name 


Center Wavelength 


a 


flux 




A 


km • 


10- erg • s-^ • cm-2 


[NeV]A3426A 


3426.7 


125.3±4.3 


118.4±3.8 


[OII]A3727A 


3728.8 


175.6±2.8 


261.9±4.1 


fNpTTTl ASSfiQA 


3SfiQ Q 


114.2±2.9 


1 1 7 Q-|_2 Q 


[NeIII]A3967A 


3969.4 


141.3±9.8 


39.8it2.6 


H<5 


4103.2 


103.3±10.6 


24.2±2.3 


H7 


4341.9 


105.1±5.9 


53.4±2.9 


[OIII]A4363A 


4364.9 


116.4±20.9 


15.8±2.7 


HeIIA4686A 


4687.6 


117.9±6.4 


54.6±2.9 


H/3 (narrow) 


4863.4 


105.5±2.1 


139.8±2.5 


H/3 (broad) 






/(Ha:i,road)/6.52 


[OIII]A5007A (normal) 


5008.9 


106.6±0.4 


1402.6±11.5 


[OIII]A5007A (broad) 


5008.8 


282.3±17.1 


150.6±11.1 


[NII]A6583A 


6586.3 


119.6±1.8 


289.7±4.3 


Hq (narrow) 


6565.6 


112.5±0.8 


540.8±3.8 


Hq (broad) 


6563.0 


2655.2±22.1 


2177.8±27.9 


[SII]A6716A 


6719.3 


122.5±3.3 


119.7±2.8 


[SII]A673lA 


6733.8 


129.6±4.2 


100.4±2.7 


(luiil)l('-gaussiaii coiiiixjiKuils iur l)r()a([ lid 




6514.6 


1283.8±69.5 


508.3±79.7 




6590.8 


2027.2±84.6 


1528.2±91.8 


(T for each line is from the 


gaussian fitting. 







Line parameters of broad H/3 arc estimated by the broad component of Ha. 
f{Ha) represents the flux of broad Ha. 

The ijaiariieters of broad component of Ha axe from fltted results by one broad gaussian 
function for broad Ha. 

The parameters of [NH]A6583A are from fitted results by two gaussian functions for broad 
Ha. 

There are two components of [OIII] A4959, 5007A doublet, one normal with not much different 
line width as that of other narrow emission lines, one broad component with broader line 
width than that of other emission lines. 



Table 2. Radio flux of SDSS J1130+0058 



Name 






ra 




dec 


int flux at 20cm 


ofl'set 
















mjy 


arcminute 


FIRST J113021.6+005820 


11 


30 


21.6330 


+00 


58 


20.000 


152.94 


0.075 


FIRST J113020.9+005828 


11 


30 


20.9480 


+00 


58 


28.210 


292.24 


0.144 


FIRST Jl 13022.4+005820 


11 


30 


22.4490 


+00 


58 


20.740 


77.58 


0.262 


FIRST J113020.1+005827 


11 


30 


20.1810 


+00 


58 


27.110 


43.96 


0.315 
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Figure 1. The observed spectrum of SDSS J1130-I-0058. 




Wavelength (A) 



Figure 2. The spectrum of SDSS J1130-f0058. The upper one 
is the observed spectrum. The middle two are the components of 
stellar lights and the power continuum. The bottom one is the line 
spectrum after the subtraction of stellar lights and continuum. 
The detailed results for CaIIA3934, 3974A are shown in top right 
corner, the dotted line represents the observed spectrum, solid 
lines is the best fitted results and line spectrum in this region. 
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Figure 3. The best fitted results for Ha. Thin solid line repre- 
sents tlic line spectrum. Tick solid line is the best fitted results. 
The dashed lines arc the two broad components for broad Ha. 
Dot-dashed lines are the narrow emission lines. The fitted broad 
Ha is shown as solid line under the observed line spectrum. 
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Figure 4. The best fitted results for H/3. Dotted line represents 
the line spectrum. Solid line is the best fitted results by scaling 
of the observed broad Ha. 
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Figure 5. The radio image at 20 cm of SDSS J1130-f0058. 
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Figure 6. The fitted results for stellar features near 
MgIbA5175A. The dotted line represents the observed spectrum. 
The solid line represents the best fitted results with ~ 0.97 by 
eight eigenspectra from spectra of stars. 
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Figure 7. The best fitted results for double-peaked broad Ha 

by elliptical accretion disk model. The dotted line represents the 
line spectrum, solid line is the best fitted results for broad Ha. 
The lower pamiel shows the line spectrum after the subtra<;tion 
of broad component of Ha. 



